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We study the largest Lyapunov exponent A and the finite size effects of a system of N fully- 
coupled classical particles, which shows a second order phase transition. Slightly below the critical 

■ energy density Uc, A shows a peak which persists for very large A'^- values (A'' = 20000). We show, 
I both numerically and analytically, that chaoticity is strongly related to kinetic energy fluctuations. 

■ In the limit of small energy, A goes to zero with a Af-independent power law: A ^ \/f7. In the 
fSJ ' continuum limit the system is integrable in the whole high temperature phase. More precisely, the 

behavior A ~ N~^^^ is found numerically for U > Uc and justified on the basis of a random matrix 
approximation . 

^ ; PACS numbers: 05.70.Fh,05.45.+b 

Recently, the interest in phase transitions occurring in finite-size systems and the study of the related dynamical 
features has stimulated the investigation of the so far obscure relation between macroscopic thermodynamical proper- 
' ties and microscopic dynamical ones. In this respect several papers appeared in the recent literature in various fields 
C — , ranging from solid state physics to lattice field theory and nuclear physics PJTc|], where there is presently a 
' lively debate on multifragmentation phase transition |p|-p^ . The general expectation is that there is a close connection 
(~| , between the increase of fluctuations at a phase transition and a rapid increase of chaoticity at the microscopic level. In 
several pioneering papers a different behavior of the Largest Lyapunov Exponent (LLE) A was found, according to the 
' O ^ order of the transition In particular, a well pronounced peak in LLE has been found for second order phase 

Q transitions, while a sharp increase has been seen for first order phase transitions. In the former case some universal 
d features have also been found, i.e. different systems show the same behavior when properly scaled pO[ |. In order to 
connect dynamical properties of systems of size N to bulk phase transitions one has to explore the continuum limit, 
N — > oo. This unfortunately is not always possible due to computer time limitations, and has been done very rarely. 
In this letter we present numerical investigations of the A^-dependence of the LLE up to N=20000, a size for which 
we already observe a certain convergence to the continuum limit. We have investigated a toy model consisting of N 
5J] classical particles moving on the unit circle and interacting via long-range forces |^ . This model shows a second order 
5^ , phase transition from a clustered phase to a homogeneous one at Uc = {E/N)c = 0.75 [||. Some results for the LLE 
of systems of moderate sizes (A^ ~ 100) have already been published [gj. The model, though relatively simple, has 
very general properties which enable us to explore the connections between phase transitions and dynamical features 
in finite systems. In particular, it could be relevant for nuclear multifragmentation where one has 100-200 particles 
interacting via long-range (nuclear and Coulomb) forces [ pd| . In this latter case a very similar caloric curve has been 
observed [ p^ and critical exponents have been measured experimentally [l^ ]. 
The main results of this letter are: 

1) The system is strongly chaotic just below the canonical transition energy Uc- The peak in A([/), found in Q for 
small systems, persists as A^ — > cxd. 

2) The increase of the LLE is related to the increase of kinetic energy ffuctuations. 

3) For U ^ 0, A ^ as U" where the exponent is found to be a = 0.5. Essentially no dependence on the system size 
is observed in this regime. A similar result was found for other systems [ p^jp^ . 

4) For [/ > J7c, A ^ as N~3. This behavior is explained by means of a random matrix approximation jist . 

5) Long-living quasi-stationary states are found in the critical region. These states look very similar to those recently 
obtained in ]l^ and simulate a discrepancy between the canonical and the microcanonical ensemble very similar to 
that one found in refs. ||l^ and more recently by other authors [p^|-^. The fact that they appear near a second order 
phase transition might be related to critical slowing down. 
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The Hamiltonian we consider is the following 



H{q,p) = K + V, (1) 



where 

K 



N N 



2"' 27V 

i—l — 1 



are the kinetic and potential energies. The model describes the motion of N particles on the unit circle: each particle 
interacts with all the others. One can define a spin vector associated with each particle uii = (cos(gi), sin{qi)) . The 
Hamiltonian then describes N classical spins similarly to the XY model, and a ferromagnetic or an antiferromagnetic 
behavior according to the positive or negative sign of e respectively In the following we will consider only 

the ferromagnetic (attractive) case and in particular e — 1. Results concerning the case e — — 1 will be discussed 
elsewhere The order parameter is the magnetization M, defined as M = -i- ™i — {Mx, My) . It is convenient 

to rewrite the potential energy V as 

y = ^(1 - (Af2 + M^)) = y (1 - M') . (3) 

The equations of motion can then be written as 

^qi=Pi , ^Pi^ ~sin{qi)Mx + cos{qi)My . (4) 



dV ^ ' dt 
In order to calculate the LLE one must consider the limit 



A = lim -In^ (5) 
t^oo t d{Q) ^ ' 



with d(£) = \l^^^i{5qi)'^ + (^Pj)^ the metric distance calculated from the infinitesimal displacements at time t. 
Therefore, one has to integrate along the reference orbit the linearized equations of motion 

. . d . ^ d'^V 



-Sq. = 5p. , -5p.. = - ^ , (6) 



J 

where the diagonal and off-diagonal terms are 



cos{qi)Mx + sin{qi)My - — (7) 



=-^-^^-9.) , . (8) 

Expression (|^) can also be written for convenience as: 

d'^V 1 

— =Mcos{q,--^)-- (9) 

where $ is the phase of M . We have integrated Eqs. (Q), (H) using fourth order symplectic algorithms with a 
time step At = 0.2, adjusted to keep the error in energy conservation below = 10^^. The LLE was calculated by 
the standard method of Benettin et al. . The average number of time steps in order to get a good convergence 
was of the order 10^. We discuss in the following numerical results for system sizes in between N=100 and N=20000. 

In fig. 1 we plot the caloric curve, i.e. the temperature as a function of U, and we compare it with the theoretical 
canonical prediction |^ (in the inset we show the magnetization). Simulations performed starting from equilibrated 
initial data, which are Gaussian in momenta at the given canonical temperature, agree very well with canonical 
predictions. In fact, it is possible to solve the stationary Vlasov equation, which represents the system in the N oo 
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limit, and obtain, under the factorization hypothesis for the probabiUty distribution, P{q,p) = f{q)g{p), and assuming 
g{p) to be Gaussian 

1 Mcos(g-$) 
= 2./o(M/T)^"P( T ) 

In the latter, Iq is the modified Bessel function of zero order and M the canonical equilibrium magnetization. The 
equilibrium probability distributions found numerically are in fair agreement with these theoretical predictions. How- 
ever, around the critical energy, relaxation to equilibrium depends in a very sensitive way on the initial conditions 
adopted. When starting with "water bag" initial conditions, i.e. a flat probability distribution of finite width centered 
around zero for g{p), and putting all particle positions qi at zero, we find quasi-stationary (long living) nonequilibrium 
states. These states have a lifetime which increases with N, and are therefore stationary in the continuum limit. We 
plot in fig. 1 the caloric curve for these states in the case N=20000. The points plotted are the result of an integration 
of 0.5 • 10^ time steps. They are far from the equilibrium caloric curve around Uc, showing a region of negative 
specific heat and a continuation of the high temperature phase (linear T vs. U relation) into the low temperature 
one. It is very intriguing that this out-of-equilibrium quasi-stationary states indicate a caloric curve very similar to 
that one found for first order phase transitions in Refs. |p7|-pl|. In that case, however, the corresponding states are 
stationary also at finite N. The coexistence of different states in the continuum limit near the critical region is a 
purely microcanonical effect. It arises after the inversion of the t oo limit with the TV — > oo one and could be 
considered as the typical signature of critical slowing down. 

We have studied how finite-size effects influence the behavior of the LLE. In fig. 2(a) we plot A as a function of U for 
various N values. In the limit of very small and very large energies, the system is quasi-integrable, the Hamiltonian 
reducing to that of weakly coupled harmonic oscillators in the former case and to that of free rotators in the latter. In 
the region of weak chaos, for U < 0.25, the curve has a weak A^-dependence. Then A changes abruptly and a region 
of strong chaos begins. In Ref. § it was observed that in between U = 0.2 and U = 0.3 a different dynamical regime 
sets in and particles start to evaporate from the main cluster. A similar regime was found in Ref. [lO| ]. This behavior 
is also similar to that found in Ref. at the solid-liquid transition. In this region of strong chaoticity we observe a 
pronounced peak already for N — 100 1^, which persist and becomes broader for N ~ 20000. The location of the 
peak is just below the critical energy at [/ ~ 0.67 and depends very weakly on N. At variance with what suggested in 
Ref. [|o| the peak does not grow with N . 

The standard deviation of the kinetic energy per particle a{K)/^/N is plotted in fig. 2(b). In the low energy region 
this quantity is in agreement with the canonical calculation (j{K) / ^/N oc U . In correspondence to the Lyapunov peak, 
we observe also a sharp maximum of kinetic energy fluctuations, though finite size effects are stronger for the LLE 
than for kinetic energy fluctuations. Thus, probes of chaotic behavior (LLE) and thermodynamical quantities (e.g. 
kinetic energy fluctuations) seem to be strongly related. We discuss here an intuitive interpretation of the relation 
between LLE and kinetic fluctuations. Each of the linearized equations (^ contains diagonal (|^) and off-diagonal (||) 
terms. Since the off-diagonal terms result from a sum of incoherent terms, we can, in a first approximation, neglect 
them. The diagonal term is of order (see Eqs. (^,(||)) and averages to < >= T + 1 — 2U. If this term 
would be constant in time, the LLE would be zero. In fact in this case one gets the equations for uncoupled harmonic 
oscillators. However, there are fluctuations, which give a non-zero coupling, whose standard deviation a{M'^) is related 
to the one of the kinetic energy, a{M'^) = 2a{K)/N, considering the relationship T = 2 < K > /N. This indicates 
that the LLE is strictly related to kinetic energy fluctuations, but this relation is not simple and quite difficult to 
extract analytically (some indications in this sense were recently proposed also in Ref. ) . In the low energy phase 
{U < 0.25) it is possible to work out a more stringent relation. In this case, the components of the tangent vector 
sum up incoherently to give for the average growth a term of the size ^/NIVP. It is then quite natural to associate 
the Lyapunov exponent to the inverse time-scale given by the fluctuations of the average growth 

A^ ^ a(VTVA/2) ^ 2^ , (10) 



Then, substituting the canonical estimation for kinetic energy fluctuations in Eq. (p_0D, we get A oc vJ7. We have 
tested numerically this prediction (see flg.3(a)). The 1/2 power law at small [/-values is fully confirmed. We have 
checked numerically that off-diagonal terms (^ cannot be completely neglected - expecially in the strong chaotic 
region. This latter important remark is also relevant for the application of a recently derived formula for the LLE 0| 
(see also [p2lp5|). A similar power law behavior was also found for other systems p6[. 
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At variance with the N-independent behavior observed at small energy U, strong finite size effects are present above 
Uc- In fig. 3(b) we show, for U > Uc, how the LLE goes to zero as a function of N. We also plot in the same figure 
a calculation of the LLE using a random distribution of particle positions qi on the circle in the equations for the 
tangent vector (|^). The agreement between the deterministic estimate and this random matrix calculation is very 
good. We find also that A scales as 7V~3 , as indicated by the fit in fig. 3(b). This can be explained by means of an 
analytical result obtained for the LLE of product of random matrices . If the elements of the symplectic random 
matrix have zero mean, the LLE scales with the power 2/3 of the perturbation. In our case, the latter condition is 
satisfied and the perturbation is the magnetization M. Since AI scales as N~^, we get the right scaling of A with 
N. This proves that the system is integrable for J7 > C/c as iV ^ oo. This result is also confirmed by a recent more 
sophisticated theoretical calculation . 

In conclusion, we have investigated the Lyapunov instability for a system with long-range forces showing a second 
order phase transition. We found strong finite size effects in the LLE. The LLE is peaked just below the critical 
energy, where kinetic fluctuations are maximal. Away from the transition region, the LLE goes to zero with universal 
scaling laws which can be explained by simple theoretical arguments. We think that this toy model contains all the 
main ingredients to understand the general behavior of the LLE in more realistic situations. 
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FIGURE CAPTIONS 

Fig.l Theoretical predictions in the canonical ensemble (full curve) for T vs. U in comparison with numerical 
simulations (microcanonical ensemble) for N=100,1000, 5000, 20000. The vertical line indicates the canonical critical 
energy Uc = 0.75. We plot also the microcanonical results for the Quasi-Stationary States (QSS) in the case N=20000 
(full circles). In the inset we show the magnetization vs. U, again the full line is the canonical theoretical prediction. 

Fig. 2 (a) Numerical calculation of the largest Lyapunov exponent as a function of U for various system sizes: 
N=100,1000, 5000 and 20000. (b) Kinetic energy fluctuations vs. U. The vertical line indicates the canonical critical 
energy Uc = 0.75 

Fig. 3 Behavior of the largest Lyapunov exponent (LLE) for U much smaller (a) and much greater (b) than 
Uc — 0-75. In panel (a) LLE shows a universal law which can be fitted by a 1/2 power law (full line). No N- 
dependence is found. In (b) the LLE - for different N and energies - is compared with a calculation done with a 
random choice of particle positions (diamonds). The latter follow a power law with an exponent -1/3 (dashed line) 
(see text). 
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Fig. 1 revised 

Latora, Rapisarda & Ruffo - Lyapunov instability and finite size... 
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Fig. 2 revised 

Latora, Rapisarda & Ruffo - Lyapunov instability and finite size. 
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Fig. 3 revised 

Latora, Rapisarda & Ruffo - Lyapunov instability and finite size. 
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